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Abstract 

Is it allowed, in the context of the Lagrange multiplier formalism, to assume that nonholo- 
nomic constraints are already in effect while setting up Lagrange's function? This procedure 
is successfully applied in a recent book [L. N. Hand and J. D. Finch, Analytical Mechanics] 
to the problem of the rolling penny, but it does not work in general, as we show by means 
of a counterexample. It turns out that in many cases the use of nonholonomic constraints in 
the process of construction of the Lagrangian is allowed, but then the correct equations of 
motion are the little known Voronec's equations. 



1 



In the Lagrange multiplier formalism as applied to nonholonomic systems, the Lagrangian is 
written as if there were no constraints. The nonholonomic constraints are taken into account in the 
formulation of the equations of motion, but not during the construction of the Lagrangian. Setting 
up the Lagrangian assuming that the constraints are already in effect is completely equivalent to 
substituting the constraint equations into the Lagrangian written as if there were no constraints. 
It is tempting to take it for granted that the ensuing reduced Lagrangian together with the relevant 
constraint equations always lead to the correct equations of motion for the system. The procedure 
just described is successfully used in a recent book [1] to solve the problem of a penny rolling on 
an inclined plane. Unfortunately, contrary to what the mentioned book appears to suggest, this 
approach is not valid in general, as we proceed to show with the help of a counterexample. 

Consider a homogeneous sphere rolling without slipping on a horizontal plane. This problem 
is treated by the Lagrange multiplier method in [2]. Let X, Y, Z be cartesian axes fixed in space 
with the Z-axis perpendicular to the plane. The principal moments of inertia with respect to the 
center of the sphere are all equal to 2mR 2 /5. With x, y the coordinates of the center of the sphere, 
the Lagrangian, being equal to the kinetic energy, is given by 



L = ^ + f) + ^f^ . (1) 



The constraint equations are 



x = Ru)y = R9 sin <f> — Rip sin 9 cos , y — —Ruj x = —R9 cos <p — Rip sin 9 sin <p . (2) 
In terms of the Euler angles 0, 9, ip the Lagrangian (1) takes the form 

L = ^(i 2 + y 2 ) + ^!(0 2 + ^ + ^ 2 + 2^cos^) . (3) 
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According to the method employed in [1], which takes into account the rolling constraint in 
the construction of the Lagrangian, the kinetic energy is written in terms of the rotational degrees 
of freedom alone by taking the moments of inertia with respect to the contact point of the sphere 
with the plane. The Lagrangian is now 

- 1 7mR 2 . , o x 1 2mi? 2 „ . . 

L=2— + + a£ ■ (4) 

This is exactly the Lagrangian that one obtains by inserting the constraint equations (2) into the 
Lagrangian (1), which was written as if there were no constraints. In terms of the Euler angles the 



reduced Lagrangian (4) becomes 



1 = 7 -^- {9 2 + i? sin 2 9) + ^- (0 2 + i? cos 2 6 + 2U cos 6) . (5) 
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Since the variables x and y do not appear in L, according to the reasoning in [1] the constraint 
equations (2) are no longer relevant to the formulation of the equations of motion. In particular, 
the Lagrange equation for 9 is 

— ( — ^ ) - — - = 7#-5?/> 2 sm#cos# + 20?/>sm# = . (6) 

dt \d9 J o9 

The treatment of this problem by the Lagrange multiplier method shows that the two Lagrange 
multipliers vanish [2]. According to equation (58d) of [2] the correct equation of motion for 9 is 

+ 0V>sin0 = O , (7) 

which is completely different from equation (6). The differential equations (6) and (7) generally 
yield different solutions for 9 because <f>, 9, ip, 0, 9, ip can be arbitrarily chosen at any particular 
instant t . It is clear, therefore, that the approach suggested in [1] lacks generality, since it works 
for the rolling penny but fails for the rolling sphere. 

In a previous paper [3] we remarked that it is possible to perform a reduction of the Lagrangian 
taking into account the constraints, but in this case the correct equations of motion are Voronec's 
equations. Given a dynamical system described by the configuration variables q±, . . . , q n , suppose 
the first m velocities are independent and the k = n — m remaining velocities can be expressed in 
terms of the independent ones by means of the equations 

m 

q m +i - a ijQj = , l — l,...,k , (8) 

where the coefficients aij are functions of the generalized coordinates q±, . . . , q n . 

Let L be the Lagrangian written without taking into account the nonholonomic constraint 
equations (8). If the last k velocities are eliminated from the Lagrangian by means of equations 
(8), a reduced Lagrangian L results: 

L(q u . . . , q n , q u . . . , q n , t) = L(q u . . . , q n , q u . . . , q m , t) . (9) 
Voronec's equations of motion are [3, 4] 
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d (8l\ 81 * dl * ™ dL 

where 

, „ _ da vi da vj A / da ui da vj \ 

13 ~ dq 3 d qi + ^\dq m+ ^ dq m+ ^) ■ W 
Setting qi = 0, g 2 = 9, q 3 = ip,q4 = x, q$ = y, in the present case m = 3 and k = 2. The 
constraint equations (2) can be written in the form (8) with 



an = , ai2 = R sin , a 13 = —R sin cos , a 2 i = , a 2 2 = — -R cos , a 2 3 = — R sin # sin . 

(12) 

The definition (11) furnishes immediately the only nonvanishing coefficients b\-\ 



b\ 2 = — &21 = — -Rcos0 , b[ 3 = —bli = — Rsm6sm(f) , 623 = ~^32 = -Rcos#cos0 , (13) 

b\ 2 = — &21 = — -Rsin0 , 6^3 = —bli = Rsin9 cos0 , 633 = — ^32 = -R cos 6* sin . (14) 
It follows that Voronec's equation for 9 is 

ImR? • IsmR? 

— - — 9 — mR 2 ip 2 sin 9 cos 9 H — <j)ip sin6> = mx^^ip + + my(\? 21 § + &23VO • (15) 

o o 

The use of the constraint equations (2) and a little algebra reduce the above equation to 

+ 0V>sin0 = O , (16) 

which coincides with the correct equation for 9 furnished by the Lagrange multiplier method. It is 
equally straightforward to check that the remaining Voronec equations for and ip coincide with 
those obtained by the Lagrange multiplier method. 

In short, the method employed in [1] is not valid in general and should not be taught to students. 
The correct result obtained for the rolling penny is the product of a mere accident. 

As a general rule, the use of nonholonomic constraints while setting up the Lagrangian is 
allowed, but the correct equations of motion are Voronec's equations, and not the ones given by 
the Lagrange multiplier method. 
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